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Abstract. 

We find that AGB stars separate in the 25-12 vs. 12-K color-color 
diagram according to their chemistry (O, S vs. C) and variability type 
(Miras vs. SRb/Lb). While discrimination according to the chemical 
composition is not surprising, the separation of Miras from SRb/Lb vari- 
ables is unexpected. 

We show that "standard" steady-state radiatively driven models pro- 
vide excellent fits to the color distribution of Miras of all chemical types. 
However, these models are incapable of explaining the dust emission from 
O-rich SRb/Lb stars. The models can be altered to fit the data by pos- 
tulating different optical properties for silicate grains, or by assuming 
that the dust temperature at the inner envelope radius is significantly 
lower (300-400 K) than typical condensation temperatures (800-1000 K), 
a possibility which is also supported by the detailed characteristics of 
LRS data. While such lower temperatures are required only for O- and 
S-rich SRb/Lb stars, they are also consistent with the colors of C-rich 
SRb/Lb stars. 

The absence of hot dust for SRb/Lb stars can be interpreted as a 
recent (order of 100 yr) decrease in the mass-loss rate. The distribution 
of O-rich SRb/Lb stars in the 25-12 vs. K-12 color-color diagram shows 
that the mass-loss rate probably resumes again, on similar time scales. 
It cannot be ruled out that the mass-loss rate is changing periodically 
on such time scales, implying that the stars might oscillate between the 
Mira and SRb/Lb phases during their AGB evolution as proposed by 
Kerschbaum et al. (1996). Such a possibility appears to be supported by 
recent HST images of the Egg Nebula obtained by Sahai et al. (1997), 
the discovery of multiple CO winds reported by Knapp et al. (1998), and 
long-term visual light-curve changes detected for some stars by Mattei 
(1998). 



1. Introduction 

AGB stars are long-period variables (LPV) and show a variety of light curves. 
Based on visual light curves, the General Catalog of Variable Stars (GCVS) 
defines regular variables, or Miras, semiregular (SR) variables, and irregular 
variables (L). The distinctive features are the regularity of light curves, their 



1 



amplitude and their period. Some types are further subdivided based on similar 
criteria (e.g. SRa, SRb, SRc,...). 

AGB stars of different variability types cannot be distinguished by consid- 
ering IRAS PSC data alone. Kerschbaum, Hron and collaborators (Kerschbaum 
& Hron 1996, and references therein; hereafter KH) have studied IR emission 
of various types of LPVs by combining IRAS broad-band fluxes with near-IR 
observations and spectra from the IRAS LRS database. They find that SRa 
variables appear to be a mixture of two distinct types: Miras and SRb variables. 
SRb variables have somewhat higher stellar temperatures and smaller optical 
depths than Miras, and can be further divided into "blue" and "red" subtypes, 
the former showing much less evidence for dust emission than the latter. 

Particularly intriguing are differences between the sLRS spectra of SRb 
variables and Miras with "10" /xm silicate emission feature (LRS class 2n). The 
peak position of the "10" fim feature for SRb sources is shifted longwards, rela- 
tive to the peak position for Miras, by 0.2-0.3 /im. In addition, the ratio of the 
strengths of silicate emission features at 18 ^m and 10 /im, F\%/F\q, is larger 
for SRb variables than for Miras. These results are confirmed at high statistical 
significance by Marengo (1998), who determines the "10" //m feature peak posi- 
tion by fitting a polynomial to the 9-11 /jm spectral region. He finds a difference 
in median peak positions of 0.2 /im, with probability larger than 0.9999 that the 
two distributions are different. 

Motivated by these results, we further analyze correlations between IR emis- 
sion and variability type for AGB stars, and their implications for the models 
of steady-state radiatively driven outflows. We base our analysis on IRAS data 
and a set of JHKLM photometric observations for about 600 LPVs, kindly made 
available to us by Franz Kerschbaum. 

2. Distribution of AGB stars in K-12-25 Color-color Diagram 

Following the results of KH, we consider only SRb, Lb, and Mira variables, and 
further divide them according to grain chemistry. Fig. 1 displays the distribution 
of sources in the K-12-25 diagram []. 

Stars in our sample clearly separate according to grain chemistry (O vs. C 
stars, S stars cannot be distinguished from O stars). Such a clear separation 
is not very surprising since the optical properties of silicate and carbon grains, 
which by and large determine the spectral shape of dust emission, are signif- 
icantly different at these wavelengths. Stars also separate according to their 
GCVS variability type into two well defined regions, Miras and non-Miras (i.e. 
SRb and Lb, which exhibit similar behavior). This is in sharp contrast with the 
distribution of the same sample in the IRAS 12-25-60 diagram, which shows no 
such division. 

A scheme which classifies stars into 4 groups according to their variability 
type (Miras, m, vs. non-Miras, n) and grain chemistry (O stars with silicate 
grains and C stars with carbonaceous grains) is shown in Fig. 1 by solid lines. 



1 We use colors denned as A2 — Ai = \og(F v (\2) / F„(\i)) , where F v {\\) and F„{\2) are fluxes at 
wavelengths Ai and A2. 
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Figure 1. K- 12-25 color-color diagram for mass losing AGB stars. 
Various symbols correspond to source variability type and grain chem- 
istry, as marked. Solid lines display classification into 4 regions accord- 
ing to grain chemistry (O and S stars with silicate grains vs. C stars 
with carbonaceous grains), and variability type (Miras, m, vs. non- 
Miras, n). Dashed line separates On sources near the Rayleigh- Jeans 
point (RJ) which have no, or very little dust. 



The boundaries, required to be straight lines, are determined by maximizing the 
efficiency and reliability of the resulting classification. The dashed line separates 
On sources around the Rayleigh-Jeans point and is determined by assuming a 
comparable scatter in both colors. 



3. Interpretations of the Source Distribution in K-12-25 Color-color 
Diagram 

3.1. Predictions of the "Standard" Models 

We first attempt to model the source distribution by employing spectra expected 
for a steady-state radiatively driven outflow with either silicate or amorphous 
carbon grains (Ivezic & Elitzur 1995, and references therein). We define "stan- 
dard" models by 2500 K black body stellar spectral shape, and dust temperature 
at the inner envelope edge of 1200 K for carbon grains and 700 K for silicate 
grains. These two temperatures mark the "standard" model tracks in Fig. 2. 
For each type of grain, the model track starts at the Reyleigh-Jeans point, with 
the positions along the track parametrized by dust optical depth. 

The track for carbon grains passes through the distribution of stars with 
carbon dust, i.e. through the regions Cn and Cm. The scatter of points is 
marginaly consistent with the photometric errors, and non-Mira stars appear to 
be smaller optical depth counterparts of Miras. The model track for silicate dust 
passes through the distribution of Miras in Om region, although not through the 
region with the highest density of sources. The offset of about 0.5 mag (0.2 on 
plotted log scale) is somewhat large, but still not sufficient to clearly rule out 
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Figure 2. K-12-25 color-color diagram for mass losing AGB stars. 
Symbols correspond to source variability type and grain chemistry, as 
marked in Fig. 1. Long-dashed lines are model tracks for silicate 
grains with various dust temperatures at the inner envelope radius, 
Ti, as marked in the figure. Models for amorphous carbon grains are 
displayed by short-dashed lines. "Standard" models have Ti=700 K 
(silicate grains) and Ti=1200 K (carbon grains). 



the model. However, the model prediction is in dire disagreement with the 
distribution of non-Miras with silicate dust (region On). For example, for 25-12 
color of -0.3, the discrepancy in 12-K color is about 0.75 on log scale, or almost 
2 magnitudes. 

3.2. Variations of the "Standard" Models 

There are several ways to augment the "standard" models without significantly 
changing the paradigm of steady-state radiatively driven outflow. 

Stellar spectral shape can effect the model spectrum when dust optical 
depth is not large. For example, increasing the stellar temperature moves a 
model track to the left in the K-12-25 color-color diagram, as required by the 
data. However, we are unable to produce satisfactory model tracks by varying 
the stellar temperature within a plausible range (< 4000 K). Another possibility 
is that a black-body spectral shape is not an adequate description of the true 
stellar spectrum. This can be ruled out by considering the K-L-M color-color 
diagram: differences in median near-IR colors between Miras and non-Miras are 
about 3-4 times smaller than required to account for the observed differences in 
12-K color. Such small differences in near-IR data also rule out an absorption 
feature in the K band, e.g. H2O feature as seen in the data by Matsuura et al. 
(1998), as an explanation for redder 12-K colors observed for Miras. Further- 
more, similar separation of stars into 4 regions can also be seen in the L-12-25 
color-color diagram. 

Models in agreement with data for non-Miras can be produced by altering 
the spectral shape of absorption efficiency for silicate grains. There are two 
required changes: the ratio of the strengths of silicate features at 18 /jm and 
10 fim has to be increased by a factor of about 2-3, and the peak position of 
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Figure 3. Model spectra in the LRS region obtained for silicate grains 
and 2 values of the dust temperature at the inner envelope radius, Ti 
(solid and dashed lines). The dotted line shows the stellar contribution 
to the spectra. Note the horizontal shift between the features, although 
both are calculated with the same grain optical properties. 

the "10" micron feature has to be shifted longwards by about 0.2-0.3 /urn. Such 
absorption efficiency results in a model track which leaves the Rayleigh-Jeans 
point at a larger angle than before, and passes exactly through the observed 
source distribution. 

Another way to make "standard" models agree with the data is to change 
Ti, the assumed dust temperature at the inner envelope edge. The model tracks 
for (unaltered) silicate grains with Ti = 300, 500 and 700 K are shown in Fig. 
2 by solid lines. They demonstrate that the distribution of stars with silicate 
dust in the K-12-25 color-color diagram can be reproduced by varying Ti in the 
range 300-700 K, without any change in the adopted absorption efficiency for 
silicate grains. It is remarkable that these models are also capable of explaining 
the differences between Miras and non-Miras found in LRS data. Fig. 3 displays 
two model spectra obtained with visual optical depth of 0.6, and Ti = 300 K 
(dashed line) and 700 K (solid line). The dotted line shows the stellar spectrum. 
Because the hot dust has been removed in the model with Ti = 300 K, both the 
blue and red edges of the silicate 10 fj,m emission feature are shifted longwards by 
about 0.2-0.3 //m, causing the whole feature to appear shifted although the peak 
position is the same in both models. In addition, this model has a larger Fig/Fio 
flux ratio, as observed. This success and the application of Occam's razor hints 
that the absence of hot dust is the most likely explanation for differences in IR 
spectra between Mira and non-Mira stars with silicate dust. 

The distribution of Miras and non-Miras with carbon dust is fairly well 
described by the "standard" model with Ti = 1200 K, and thus the assumption 
of significantly lower Ti is not required by the data. However, such models are 
not ruled out by the data as shown by model tracks for Ti = 600, 800 and 1200 
K displayed by dashed lines in Fig. 2. 



5 



4. Disussion 



Absence of hot dust, as implied by low Ti, might be a natural consequence of 
a significant recent decrease in mass-loss rate. From the distribution of sources 
in the On region, we infer a time scale of about 100 years. What happens after 
that? The distribution of On sources ends quite sharply at its blue edge. If 
envelopes continued to freely expand, distribution would extend smoothly to 
bluer 12-K, and redder 25-12 colors. The absence of these sources argues that 
mass-loss resumes after about 100 years. This idea, that stars oscillate between 
the SRb and Mira phases, was proposed by Kerschbaum et al. (1996) based on 
similar number densities and scale heights. 

This possibility seems to be supported by a number of new observations. 
Recent HST images of the Egg Nebula obtained by Sahai et al. (1997) show 
concentric shells whose spacing corresponds to a dynamical time scale of about 
100 years. The discovery of multiple CO winds reported by Knapp et al. (1998) 
is in agreement with the hypothesis of variable mass-loss. Furthermore, they 
find that the velocity and mass-loss rate for each individual component follows 
the relation v oc M 1 / 3 found for stars with single CO winds (Young, 1955), and 
theoretically expected for radiatively driven outflows (ivezic, Knapp & Elitzur, 
1998). Finally, there is mounting direct evidence that long-term visual light 
curve changes are consistent with the transition between Mira and non-Mira 
variability types (Mattei, 1998). 

If proven correct, the mass-loss variations on time scales of about 100 years 
would significantly change our understanding of the stellar evolution on the 
AGB. What type of observations could rule out this hypothesis? The most 
obvious and direct test is mid-IR imaging. When Ti decreases by a factor of two, 
the inner envelope radius increases by about a factor of four. Such a difference 
should be easily discernible for a few dozen candidate stars with largest angular 
sizes, which are already within the reach of the Keck telescopes (Monnier, 1998). 
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